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We develop a general formalism for analyzing linear perturbations in multiple-field cosmological 
inflation based on the gauge-ready approach. Our inflationary model consists of an arbitrary number 
of scalar fields with non-minimal kinetic terms. We solve the equations for scalar- and tensor-type 
perturbations during infiation to the first order in slow-roll, and then obtain the super-horizon 
solutions for adiabatic and isocurvature perturbations after infiation. Analytic expressions for power- 
spectra and spectral indices arising from multiple-field infiation are presented. 
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Observations suggest that the early Universe underwent a period of accelarated expansion called cosmological 
inflation. In addition to providing a causal mechanism for the generation and evolution of large-scale structure 
formation, inflation also leads to elegant resolutions of a number of puzzles of the Big-Bang theory, such as the 
: isotropy, horizon and the flatness problems Q, Q, |^ . The simplest implementation of inflation is achieved by assuming 
I> ' that the matter is described by a single scalar field, the inflaton The quantum fluctuations of the scalar field 

, generated during inflation become classical after crossing the event horizon, seeding the observed density perturbations. 
■ In addition, inflation also generates metric perturbation in the tensor sector, leading to a stochastic gravitational wave 
background. Indeed, the recent high accuracy CMBR data from the WMAP satellite lij do indeed support the general 
predictions of inflationary cosmology, 
iy-^ I It was, however, realized, that with a single scalar field the inflationary scenario suffers from what is called the 
, graceful-exit problem namely, achieving sufficient inflation consistent with the observed density perturbations, 

before the Universe exits from the inflationary epoch. Kofman et al. 7] suggested the possibility of using a first 
short stage of double inflation in order to generate a large value of a scalar field required for a second, longer stage 
with a graceful exit. Linde showed that one requires at least two scalar fields to overcome the graceful exit issue, 
O without modifying Einstein gravity, and without sacrificing natural initial conditions. There are other motivations 
for incorporating multiple fields in the dynamics of cosmological infiation. For example, when constructing models of 
^ inflation inspired by particle physics theories such as low energy effective supergravity derived from superstrings, one 
• • obtains many scalar fields (see Q for a recent review). This calls for a general framework for handling cosmological 
. ^ ' perturbations in a situation where the matter sector consists of an arbitrary number of scalar fields. 
^> , Cosmological perturbations in a single field inflation has been thoroughly investigated in the pas t 0, 0, 

' following the seminal paper of Bardeen In the context of multiple-field inflation, Starobinskv |l3j| obtained an 
expression for density perturbations with an arbitrary number of scalar fields interacting between themselves through 
gravity. In a consistent treatment of cosmological perturbations with more than one field, one should consider the role 
of isocurvature, or entropy modes in addition to the adiabatic, or curvature perturbations 0,0,0,0. Indeed, it 
is quite possible for the two to be correlated, leading to distinct observational results Recently Gordon et al. 0| 
analyzed the evolution of adiabatic and isocurvature perturbations in multicomponcnt inflation, where they performd 
a local rotation in the field space to separate out the adiabatic and entropy modes, while Wands et al. |23| studied 
possible observational aspects of adiabatic and isocurvature spectra produced by two-field inflation. For a somewhat 
different approach, see Malik et al. |2]J. A method for treating density perturbations in multicomponent inflation was 
proposed in [13, but see also ^.^ol^. 

In a recent paper |26| we presented a general formalism to analyze cosmological perturbations in inflation driven 
by multicomponent scalar fields using the gauge-ready method developed by Hwang and colleagues [271 Isof . 
This approach follows from a suggestion by Bardeen |31l] . that rather than imposing a particular gauge condition 
while dealing with cosmological perturbations right from the beginning, it is often advantageous to express the 
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perturbations without specifying any gauge. Thus one has the flexibihty of adopting different gauge conditions at a 
much later stage, depending upon the nature of each problem. Moreover, it becomes easy to relate results between 
various gauge-dependent and gauge-invariant techniques. 

This paper is organized as follows. In Section |21 we present the equations describing multicomponet scalar fields 
with a non-trivial field metric and having non-minimal kinetic terms coupled to Einstein gravity. By introducing a set 
of basis vectors based on the Gram-Schmidt orthonormalization technique, we can discriminate multiple-field effects 
from single-field ones. Metric perturbations are discussed in Section|31 These are conveniently decomposed into scalar, 
vector and tensor modes. Here we introduce the gauge-ready approach to cosmological perturbations and present the 
equations governing density perturbations in the gauge-ready form, as well as in terms of gauge-invariant variables. 
Next, we introduce the slow-roll variables in Section ^ and proceed to obtain the solutions to scalar and tensor 
perturbations during inflation in the first order in slow-roll. We then proceed to study perturbations after infiation. 
In Section[31we discuss adiabatic and isocurvature perturbations. We calculate the power-spectra and spectral indices 
for adiabatic, isocurvature, correlated, and tensor modes. These should be useful in comparing theoretical predictions 
of various inflationary models with observations. We conclude in Sectional 



2. THE INFLATIONARY MODEL 



2.1. Scalar fields 



In this Section, we explain our notation and set up the basic equations needed for our analysis. We consider Einstein 
gravity coupled to an arbitrary number of real scalar fields. The Lagrangian then is 

1 „ 1, 



(2.1) 





1 „ 1 



^ ( iji^ - ^g'"'a^(/.^G9,0 - y(0) 



were R is the scalar curvature, Kq = SttG, and we set c = 1. For the scalar fields we use a vector notation, = (0"), 
with the indices a, 6, c, . . . = 1,2,3, . .. ,N labelling the A^-components in field space. Further, g = det(g^i,), and 
fj,,i>, . . . denote the spacetime indices. For repeated indices, the summation convention applies. The second quantity 
within the parentheses of Eq. 12.1|l represents a non-minimal kinetic term. Such a kinetic term appears in various 
models of high-energy physics 9J. Also F((/>) is an arbitrary scalar potential. 

The scalars (f) may be interpreted as coordinates ((/)°) on a real manifold M induced with a symmetric Riemannian 
metric G having components Gab in the field space |22l |. The field metric is chosen to be positive-definite so that 
the Hamiltonian is bounded from below. The special case of minimally-coupled fields corresponds to the situation 
Gab = Sab- From the components Gab we can define the connection coefficients F^^ in the usual manner, 

^bc = -G"'' {Gbd,c + Gcd,b — Gbc,d) ■ (2.2) 

The curvature tensor on M is introduced in terms of the tangent vectors B, C, D: 

[RiB, C)Dr ^ B'C^ ^ (F?,,, - F^,,, + Fg,F^, - Fg.F^,) B^C^ D". (2.3) 
For any two vectors A and B, we define the inner product and the norm as 

A B = A^B = A^GB = A^'GabB^, 



\A\ EE ^(aTa), (2.4) 

respectively. Here is the cotangent vector such that {A^)a = A^Gba- We also introduce the covariant derivative 
Va on M acting upon a vector A as 

VbA-^ = A% + F?,A^ (2.5) 
while, the covariant derivative on A with respect to the spacetime x'^ is 



(2.6) 
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It should be noted that the covariant derivative reduces to the ordinary derivative when it acts upon a scalar. 
By varying the action H2.1|l with respect to g^i, and ^, we obtain the gravitational field equation, 

^Gt =n = d^^cj, ■ d,cj> - <5^: • dx(t> + , (2.7) 

and the equation of motion for the scalar fields, 

g^'' {vX~^Udx4>-G'^V^V = Q, (2.8) 

where and are Einstein and energy-momentum tensors. 

It is often convenient to represent the scalar fields as effective fluid quantities. The energy-momentum tensor can 
be covariantly decomposed into fluid quantities using a time-like four- vector normalized as u'^u^ = — 1: 

Tap = ^^UaU|3 + phaf} + qaUf3 + QpUa + TTq^, 

TTaf3 = T^shlh^p - pKp. (2.9) 

Here /i, p, qa, and tTq^ are the energy density, pressure, energy flux, and anisotropic pressure, respectively; hap = 
dap + UaUp is a projection tensor based on Ua vector, (/qM" = = T^afj, and tt^ — 0. Thus, for a multicomponent 
scalar field, the above decomposition gives 

M=^I0P + ^, P=l\<P\'-V, 

qa=0 = TTap. (2.10) 

Equations 12.7|l and (|2.8() with H2.10|l provide the fundamental expressions required for describing cosmological infla- 
tion. 



2.2. Basis vectors 



It will now prove useful to introduce a set of basis vectors generated using Gram-Schmidt orthonormalization 
p^ls^ . From the vector (f> we can construct a set of N linearly independent vectors {4>^^\ <P^'^\ ■ ■ ■ , </''^''}, where, 

^(i) = 0, (n>2). (2.11) 

Let ei = be the first unit vector along the direction of the field velocity 0. Define the second unit vector 

62 to be along that part of the direction of the field accelaration Vtcp which is normal to ei so that ei ■ e2 — 0. 
Repetitively applying this Gram-Schmidt procedure generates a set of mutually orthonormal vectors {e,i}, spanning 
the same subspace as the vectors {(/)("^}. 

Introducing the projection operators Pn and P^, which project on e„ and on the subspace perpendicular to 
{ei, . . . , e„} respectively, we may then write a general unit vector e„ as. 



(2.12) 



where, 



P„ = e„et, P,^ = 1-^P„ Po^^l, (2.13) 

9=1 

and we define, 

pll=Pi=eie|, P^ = P^ = 1-PK (2.14) 

Note that when the denominator in Eq. H2.12|l vanishes, the corresponding basis vector does not exist. 

Since pH + P-^ = 1, we can decompose any vector A in directions parallel and perpendicular to the field velocity: 

A = All + = (Pll +P'^)A 

= ei(ei-A)+e2(e2- A). (2.15) 

When there is just one field, Ci by definition simply reduces to the normalized scalar | while 62 vanishes 

identically, and so do all other basis vectors. Thus the decomposition (|2.15|) enables us to distinguish between single- 
field contributions, where only Ci survives, from multiple-field ones. 
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3. THE PERTURBED UNIVERSE 
3.1. Metric perturbations 

The observed Universe is not perfectly homogeneous and isotropic. Assuming that the inhomogeneities are small 
enough, we can then treat the deviations by considering linear perturbations of the homogeneous and isotropic 
cosmological space-time described by the Friedmann-Robertson- Walker ( FRW ) model, 

ds^ = -a^ (1 + 2A) drf - 2a^B,dx' 

W (^g'lf +2C,j^dx'dx^, (3.1) 

where a{t) is the scale factor, dt = adrj, and indices . . ., run from 1 to 3 labelling the spatial components. The 

perturbed order variables A{t,x), Bi{t,x), and Cij{t,x) are based on the metric glj of the 3-surfaces of constant 
curvature K = 0,±1. Here t and ry are the comoving and conformal times respectively. We denote a derivative with 
respect to comoving time by ' = dt and one with respect to conformal time by / = 9^. The Hubble parameters in 
terms of comoving and conformal times are defined as H ~ d/a and Ti = a' /a = aH. 
Similar to the metric decomposition Eq. (|3.1|l . we can decompose the scalar field as 

4>{t,x) = (t>{t) + S<f>{t,x), (3.2) 

where the perturbation dip = {^4'°') is a tangent vector on A4, while the energy-momentum tensor is decomposed as 

T° = - fe + {^J. + p)ui] = {n+p)vi, 
a 

T)=p5]+^)^{p + 5p)S]+^% (3.3) 

The barred entities denote background variables. For notational simplicity we shall ignore the overbars unless required. 
In Eq. H3.3|l . Vi is the frame-independent flux variable, and Vi, tt^^'^^j are based on g^^"*. 
From Eqs. H2.7|l and H2.8|l . the equations for the background can be written as 



(3.4) 

H = -Uu^^ + p) + § = -l4\<Pf + ^, (3.5) 

R = 6(^2H^ + H + , (3.6) 

Vt^ + 3H(p + G-^V^V = 0, (3.7) 
/i-f 3H(/i + p) = 0. (3.8) 



We shall ignore the cosmological constant A in our work; nevertheless it can be easily included by making the 
replacements ^ ^ ^ + A/kq and p ^ p — K/ Kq. Note that we have explicitly retained K{= 0, ±1), and only at a later 
stage shall we set K — Q. 

3.2. Scalar, vector and tensor decompositions 

To make further progress, we decompose the perturbed order variables into scalar-, vector-, and tensor-type per- 
turbations. To the linear order, they decouple from one another and evolve independently. Accordingly, the metric 
perurbation variables A{t,x), Bi{t,x), and Cij{t,x) may be decomposed as 

A = a, B, = (3, + bI^\ 

a,^ggV + 7,,|,+C;:;;)+C^). (3.9) 
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The superscripts (s), {v) and {t) indicate the scalar-, vector- and tensor-type perturbed order variables. The vertical 

(3) 

bar represents a covariant derivative with respect to gl^ and the round brackets in the subscript imply symmetriza- 
tion of the indices. The scalar metric perturbations are then given by a, /3, 7 and ip. The transverse-type vector 
perturbations i?!^"* and C^"'' satisfy B^""* |j = = C''"''j. while the tensor-type perturbation C^j' is transverse-traceless 

(C''*-'*j = Both the vector and tensor perturbed order variables are based on g'i^\ We define A as a 

comoving three-space Laplacian, and introduce the following combinations of the metric variables, 

X = a(/3 + a7), k ee 3{Ha ~ ip - ^x), 

= B^"^ + aC^^'l (3.10) 

It is convenient to separate the temporal and spatial aspects of the perturbed order variables by expanding them in 
terms of harmonic eigenfunctions Q^*-"'*) (fc; x) of the generalized Helmholtz equation [TlLfl^. with k the wave vector in 
Fourier space and k = \k\. We can then write the scalar-type perturbed order variables as a(t, x) = a{t, fc)Q^*-'(A;; x), 
with similar expressions for /3, 7 and ip. The vector- and tensor-type perturbations are expanded as sj"'' = B^^^Q'f\ 

C^"'' = C*^"-* qI"', and C^j"* = C^*-'q|*'. In each of these harmonic expansions, a summation over the modes of the 
eigenfunctions is implied. Thus, the perturbed scalar fields have the expansion 

Scf){t,x) = 6ct}{t,k)Q^''\k;x). (3.11) 
Similarly, the fluid variables Vi and tt^^''^ can be expanded in terms of the harmonics as 

^(3)^ ^ ni^)Q^^) + ^^)q!^^] + ^(*)q(*);; (3.12) 
while the energy-momentum tensor in Eq. 13.3|l has the expansion 

Tl = -^j,= -{Ji + (5pi), 



T) = {p + 5p)5)+ 7r(^) q(^'; + 7r(") Q^"^ + ^(*) Q^'^] . 



(3.13) 



For a Universe having the matter sector composed exclusively of scalar fields, the quantity tt''^'' ^ in Eq. (|3.3fl vanishes 
identically. We then have to the perturbed order, 

5n = ^-Vt5(i)-a\4)\^ + VV -54), (3.14) 
5p^^-Vt5(p-a\4)\^ -VV ■5(p, (3.15) 

(/i + p)f^ = 0-(50, (3.16) 

where we have written v = w*-*-* for simplcity. It is also convenient to decompose 5p into an adiabatic part c^J/i, and 
an entropy perturbation e: 

5p = cl5n + e, (3.17) 
where = p/ fi may be interpreted as an effective sound velocity. We shall also use the notation w = pj [l. 



3.3. Gauge-Ready formalism 

We now briefly summarize the gauge-ready approach discussed in l27l I2M l29l l30l | . Under a gauge transformation, 
or coordinate shift, 5:^ — -f with (^'^, ^') = (a^^^*, a^^^'' +1^'^ ■")*), and C'""'*!* = 0, the metric and matter variables 
transform to linear order as 
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a 

A\ . 1 



\ a'^ J a 
5fL = Sfi — /i^*, 5p = 6p — p^*, S(f) — S<p — 0^*, 

It is immediately obvious from Eq. (|3.18(l that the tensor- type perturbations are gauge-invariant. For the special 
case of scalar-type perturbations to the linear order, fixing the temporal part ^* of the gauge transformation leads to 
different gauge conditions: a = (synchronous gauge), x = (zero-shear gauge), v/k = (comoving gauge), cp = 
(uniform-curvature gauge), and so on. Except for the synchronous gauge, the temporal gauge mode in the other 
gauges are completely fixed. Consequently, there is a unique correspondence between a variable in a gauge condition, 
and a gauge-invariant combination of the variable concerned and the variable used in the gauge condition. Using the 
perturbed order variables together with the variables used in the gauge condition, one can systematically construct 
various gauge- invariant variables, for example, 

k 

ip^ = ip- Hx, a^ = a-x, = v - -x, 
= dfi- fix, hx =6p- PX, 

aH a 
(fiv = •P> -j^v, dfiy = dfi~ -fiv. (3.19) 

Thus, in the zero-shear gauge, also known as the longitudinal, or conformal Newtonian gauge, we have from 
Eq. (|3.19(l . ip^ = ip, = a, and Scj)^ = Scj). Similarly, in the uniform-curvature gauge, it follows that Scj)^ = Scj) 
which in turn is equivalent to —{(p/ H)ips(p in the uniform-field gauge. In the notation of [inl | , our and ip^ correspond 
to their $ and — \E' respectively. 

Now, as is well known in the theory of cosmological perturbations, a judicious choice of gauge conditions often 
simplifies the mathematical structure of a particular problem. For example, density perturbations with hydrody- 
namical fluids are most conveniently treated using the comoving gauge, while gravitational potential and velocity 
perturbations are best handled in the zero-shear gauge. In the same spirit, the uniform-curvature gauge simplifies the 
analysis of perturbations due to minimally coupled scalar fields. Since, in general, we do not know the optimal gauge 
condition beforehand, it becomes advantageous to express the perturbations without imposing a specific temporal 
gauge condition. In other words, we write the governing equations in the gauge-ready form, which would give us the 
freedom to choose different gauge conditions, as adapted to the problem, at a later stage in the calculations. The 
equations are spatially gauge- invariant, but the temporal gauge condition remains unspecified. Once the temporal 
gauge mode is completely fixed so that no further gauge degrees of freedom are left, the resulting variables would 
then be gauge- invariant. Moreover, when a solution in a particular gauge is known, we can then easily derive the 
corresponding solution in other gauges, as well as in gauge-invariant forms. This is the basic concept of the gauge- 
ready method. The method is most useful when one considers relativistic hydrodynamic perturbations with mutually 
interacting imperfect fluids as well as kinetic components. The gauge-ready method then not only simplifies the 
analysis by enabling us to choose different gauge conditions for different aspects of the system on the fly, but also 
allows us to check the numerical accuracy by comparing solutions in different gauges. 

To implement this gauge-ready strategy, it is most convenient to derive the perturbed set of equations from the 
(3+1) ADM (33, and the (1+3) covariant [33| formulations of Einstein gravity. A complete set of these equations 
may be found in the Appendix of Ref. . In this Section we write the equations for scalar-type perturbations in the 
gauge-ready form: 



1 1 fc2 
if ^ Ha- -K + -— X- 
3 S a'' 



3K 



Hk 



(3.20) 
(3.21) 
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k^-3K So. ^a 



X + Hx - a - tp ^ Kq—t: 



is) 



k + 2iJK 



+ mvt - 



d — 



A 



X ) - 2aG-iV^y. 



fc 

+ ZH[6^ + Sp) = (/i + p) ( K — 3i?a w 



a4(/x + p) 



^ 2 fc2 - 3X 
/I + p V 3 k'^ 



-is) 



(3.22) 
(3.23) 
(3.24) 

(3.25) 

(3.26) 
(3.27) 



Equations H3.20|l - H3.27|) are the definition of k, the ADM energy constraint (Gq component of the field equation), 
the ADM momentum constraint (G^ component), the ADM propagation(GJ — ^SjG'l component), the Raychaudhuri 
equation (G* — Gg component), the equation of motion for scalar fields, energy conservation, and the momentum 
conservation, respectively. Here J/z and Sp are given by Eqs. H3.14|) and (|3.15|) respectively, while 



(3.28) 



Note that these equations are valid for any K, and for a scalar field, tt^''^ = 0. 

Equations H3.20|l - I|3.27(l . together with the background equations H3.4|l - H3.8(l . and the perturbed order variables for 
the scalar fields (|3.14|I - H3.16|I . provide a complete set of equations for analyzing scalar- type cosmological perturbations 
with multicomponent scalar fields. As we have not chosen a specific gauge so far, Eqs. H3.20|l - H3.27|l are therefore in 
the gauge-ready form. This allows us to impose any one of the available temporal gauge conditions, which would then 
fix the temporal gauge mode completely, leading to gauge-invariant variables. 



3.4. Gauge-Invariant perturbation equations 



As an illustration of the gauge-ready method, we derive some useful expressions using the gauge-invariant variables 
of Eq. (|3.19() introduced in Section lTSl These may be obtained by making judicious combinations of Eqs. H3.20|l - H3.27|l . 
Thus, from Eqs. and we obtain 



k^-3K 1 



(3.29) 



Eq. 13.23|l can be re-expressed as 



2 ^ / 



(3.30) 



Eqs. (jX^ . and lead to 



1 2/ 



fx - Ha^ = -o'«o(Ai+P)T:"x- 



(3.31) 



Eqs. (jT^ . with (jT^ yield 



6 fly + iHS^T] 

k^-3K 



(^+p)«,^ + 2i7^.(^) 



(3.32) 



8 



while Eqs. and (|X77|) 



give 



2P_ 
3 



3K TT 



P+P 



(3.33) 



Combining Eqs. (|3.29() - (|3.33() we can derive 



+ (4 + 3cl)Hip^ ~ cl^Vx 



ificl -p)- 2(1 + 3c2) 



K 



1 



3 



1 / 2ir2 



H \n + pk'^ 



(3.34) 



where we used Eq. (|3.17|l . From Eq. 



since tt'^^ = 0. Using this result the equation of motion for scalar fields becomes 



wc can draw the important conclusion that, for scalar-fields, — ~^Px^ 



2HV„ - A 



+2aVxG"^V^F, 



while Eqs. l|OT)) and become 



1 



• S(f>x 



1/ 



-x ' ' "^x — 2 
■ mtp'^ - Aipx + 2[H' + 2(7^2 _ K)] 



(3.35) 

(3.36) 
(3.37) 



where we used Eq. (|3.16|l . and the relations 



e = Sp- clSii = 5px - clSfi^, 

(1 - cl)Snx -e = 5nx- 5px = W ■ S<f>x 



(3.38) 



Eq. (|3.36|) is often called the constraint equation. These equations contain most of the physics related to inflationary 
cosmological perturbations. They are expressed in terms of gauge-invariant forms of the variables, and from Section 
13.31 we see that they retain the same algebraic forms in the zero-shear gauge. 

We note that Eq. (|3.37|l may be recast in a different way. According to Eq. H2.15(l . dip^ may be decomposed into 

components parallel and perpendicular to the field velocity, Scp^ = <5(/)| + S(p^. Using the background equation (|3.7|l . 
the constraint equation (|3.3t)|) . and the fact that — (2^))0') ' 4>\ we can write Eq. H3.37|l as 



2[n- 



101 



n' -n 



I0'l 



- 2K 



^x - ^fx 



(3.39) 



Following our discussion in Section 12.21 we know that the perpendicular component of field perturbation vanishes 
when there is only one field. In this case, the right hand side of Eq. H3.39|l vanishes, and the resulting equation is well 
known in the theory of single field inflationary perturbations |lC| . 
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4. PERTURBATIONS DURING INFLATION 
4.1. Slow-roll variables 

Continuing with our analysis, we now assume that the Universe has undergone inflation to complete flatness, so 
that henceforth we can set K = 0. This allows us to introduce the a set of functions known as the slow-roll variables: 

= r](ct>) = (4.1) 

It is also convenient to decompose r) into parallel and perpendicular components using Eq. I|2.15|l . 

The standard slow-roll assumptions are 

e = 0(C), ?/"=0(C), ^^ = 0(0, (4.3) 

for some small parameter C, with e, y^?/" and a/e?/^ much smaller than unity. If in an expansion in slow-roll variables 
we neglect terms of order O(C^), we claim that expansion to be of first order in slow-roll. Thus terms with e^, er]\ 
etc. are of second order. We list some useful relations involving the slow-roll variables: 

Vrjct)' = n\(t>'\{r, + ei) = K-'V2n^V~eiv + ei), 

nh^^Kl\<t>'\', e'^2He{e + ^^. (4.4) 

4.2. Analysis using gauge-invariant variables 

In order to solve the system of perturbation equations H3.35|l . H3.36|l and (|3.39|l . we shall find it convenient to 
introduce the variables, 

q = a (^Sc/)^ - Ji^'Px^ = « (^^(t> - ^'P^ , (4.5) 

Here q is a gauge- invariant quantity, and is a natural generalization of the single field Sasaki-Mukhanov variable [lOj | . 
We first express the constraint equation l|3.36|l in terms of the slow-roll variables 1)4. l|l as 

<^;^+H(l + e)(Px = -i«2^'.l (4.7) 

The equation for the scalar field perturbations follows from H3.35|l . 

_ (A - n^n)q = 0, (4.8) 

where 



n 



and we also used (|4.7|) . The corresponding Lagrangian follows from Eq. (|4.8|) 



(2 - e)l - 2e ((3 + e)pH + eirjt + r/ej) , (4.9) 



= \ J {'D„q^V^q + qHA-H^n)q)y/^df]d^X. 



(4.10) 
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Here i?^^^ is the determinant of the metric g^^^ of the 3-surfaces of constant curvature K = 0, see below Eq. (|3.1|l . 
The equation of motion for u is obtained by substituting its definition (|4.t)|) into Eq. (|3.39|l : 



u"-Au-—u = nri 52, 92 = 62-9, 
. H kq 1 



a|0'| V2ay^' 

Further, from Eq. 1)4. 7|l it follows that 



(4.11) 



u' + ^^u^ -qi, qi = ei-q. (4.12) 
Differentiating Eq. (|4.12|l once with respect to the conformal time and using Eq. H4.11|l . we obtain the relation 

^91) - '^^^1-2 = A^i. (4.13) 



2 V ' I/O 

We pause to note that though the equations (|4.7|l . (|4.8|) and 14.11|) have been expressed in terms of the slow-roll 
variables, they are exact, and no slow-roll approximation has yet been made. Observe that, to the leading order in 
slow-roll, the perturbation variables q and u decouple, whereas at first order, mixing between these occur. 



4.3. Quantization of density perturbations 

We now briefly discuss the quantization of the density perturbations described by the Lagrangian in Eq. 14.10|l . 
Introduce the matrix Zm„ defined as 



(^)mn = -{Z'^),nn = ^e,n ' V^E^, (4.14) 

rl 



Thus Z is antisymmetric and traceless. Expanding q — qm^m using the basis {em}, Eq. (|4.10|) may be expressed as, 

^ - i(q' + HZqfiq' + nZq) + ig^(A - -H'n)q, (4.15) 

where {il)mn — 6j,jJle„, and for notational ease, we have suppresed the indices m, n. Now redefine q using a new 
matrix R as 

q{r]) = R{'q)Q{7]), R' + TiZR^O, n^R^QR. (4.16) 

From the equation of motion H4.16|l for R, it follows that R^R and deti? are constants, so that R represents a rotation. 
Without any loss of generality, the initial value of R may be chosen as -R(?7o) = 1- Substituting the variables defined 
in Eq. (|4.16() into Eq. 14.15|l reduces the Lagrangian to the canonical form: 

^ ^^Q'^Q' + ^Q^{A-n''n)Q. (4.17) 
The corresponding Hamiltonian is then given by 

^ = in^n- ig^(A-H2fi)Q, (4.18) 

where the momentum 11 canonically conjugate to Q is 

n{T], x) = d^/dQ'^ = Q'it], x). (4.19) 

To implement the canonical quantization procedure, the variables (Q, 11) are promoted to quantum operators (Q, 11) 
satisfying the commutation relations 

[a'^Q{r],x),PQ{rj,x')] = [a'^Il{rj,x), pilir^,x')] = 0, 

[a^Qiv, x),l3Il{ri, x')] = ia^(3S{x - x'), (4.20) 
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where the delta function is normaHzed as 



6{x-x')Vg(^d^x = 1. 



(4.21) 



Here we have introduced the vectors a, /3 with components a^n, Pm in the basis {6^} to avoid writing the indices 
m, n in the commutators. Since we are considering spatiaUy flat hypersurfaces (K — 0), the operator Q may be 
expanded in a plane wave basis as 



(27r)3/2 



Qk{v)aie 



-t —ik-: 



h.c. 



(4.22) 



with a similar expansion for H. It immediately follows from Eq. (|4.1t)|) that q must now be interpreted as the operator 
q satisfying a mode expansion identical to Eq. 14.22|1 . The creation and annihilitation operators and satisfy 



[a^dk,Pak'] = [a^dl,Pal,] = 0, 
To maintain consistency of the commutation relations 14.2011 and H4.23() . the Wronskian condition 

w{Qk,Ql} = Q'MQUv) - Qk{v)Qk{v) = ^ 



(4.23) 



(4.24) 



must be satisfied. From the mode expansion (|4.22|l and the Hamiltonian (|4.18|l , it follows that the equation of motion 
for Qk is 



(4.25) 



It may be easily verified using Eq. 14.25|l that the Wronskian satisfies dW{Qk, Ql}/di] = 0. 

We also interpret the variable u introduced in Eq. (|4.6|l as an operator u, and after performing a mode expansion 
identical to that of Q in Eq. ()4.22|l . it follows from Eq. 14.11|l that the modes Uk satisfy 



Uk = Tii]^q2k, q2k = (62 • em)qk, 



(4.26) 



or, equivalently, from Eq. H4.13|l . 



1 / (1/0)' 
qik = (6i • ern)qk- 



(4.27) 



4.4. First order solution 



In order to present the solution of the scalar perturbation equations, it is convenient to introduce the time rji-i when 
the mode with wave number k crosses the Hubble radius during inflation, so that the relation 

■Him) = k (4.28) 

is satisfied for each k. Consequently, the infiationary epoch can be separated into three regions: the sub-horizon 
region (7i <C fc), the transition region (Ti. ^ k), and the super-horizon region (7i ^ fc). We now discuss each of these 
in turn. 

In the sub- horizon region, we solve Eq. 14.25|l with the Ti^il term subdominant compared to k^. The solution is 
obtained in the limit k/Ti. — > oo for fixed k as 

Q,(rj) ^ ^e^'^^^-^"\ Rim)^l. (4.29) 

Since one is usually interested in calculating quantities at the end of inflation, this region is therefore irrelevant. 
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Considering next the transition region, we introduce the time 77- when the sub-horizon epoch ends and the transition 
region begins. In a sufficiently small interval around 77-^ we can then apply slow-roll to the Eq. (|4.25|) keeeping all the 
terms, but taking the slow-roll functions to be constant to the first order. The initial conditions are chosen as 

1 i\/k 
Qfe(^-) = ^1, Q'fe('7-) = ^l, i?(r,_)-l. (4.30) 



Integrating the relation for Ti.' in Eq. (|4.4(l with the initial conditions for the transition region, we obtain 

so that T-Hji-h) — k. Differentiating 9 in Eq. 14.111) yields 9' /6 — —H{1 + e + 77II), which can be integrated with the 
result 

/ z \ (i+2^«+''«) 
9iz)^9n[—] , z = kr], 



K H-j-l 



Solving Eq. (|4.16|) for the rotation matrix R with the initial conditions H4.30|l yields 

z 



zn = k-nn- (4.32) 



R{z) = ( — ) , z_ = /c77_. (4.33) 



Since the time-dependent terms in the matrix in Eq. (|4.9|) are of first order, we can take = flijj-H) = fi-^ in the 
transition region. Then the matrix f2 is given by 

n ^ R-^{z)nnR{z) (4.34) 



with 



fin + 'S[Sn,ZH] [In— + -h-i en] , (4.35) 

z-H 4 



^('7) = -^21 " 



3 V. (1 - e) 
37^2 

Sn = Sim)- (4.36) 



el - 2 + 2e(ei • e„)(ei • e„)^ 



Here the second equality is valid to the first order in slow-roll, and we have used the notation = el^^M^en- We 
also made the assumption that those components of a?M'^ /Ti'^ which cannot be expressed in terms of the slow-roll 
variables are of first order. Because 5h and Zn arc both of first order, we can take fl = fin in Eq. H4.35() to be a first 
order quantity. 

In order to write the equation for the mode Qk in the transition region, we will find it convenient to define 
Qk = RnQk{z) and fl = Rn^Rn^, with Rn = R{zn)- From Eq. (|4.33|) . we have to the first order, Qk{z) = Qk{z), 
while from Eq. (|4.35() we conclude that fl = fin within a small region around zn ■ Using the above results in Eq. 14.25|l , 
the mode equation for Qk may be written in terms of Qk as 

Qfe, + (^a - Qfc = 0, jy^^^l + 35n. (4.37) 

This equation is similar to the one obtained for the single-field inflation, except that this is a matrix equation. The 
solution is then given in terms of the Hankel functions of matrix valued order vn, 

Qk{z) = Vi[ci(fc)i?W(^) + C2(fc)i?(^)(z)], 

i^H^l^ + Sn- (4.38) 
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We wish to match the sohition in Eq. (|4.38() so that in the hmit k/H — > oo, the modes approach plane waves, 
Qk{z) = e'^/\/2fc, see H4.29|l . For \z\ ^ 1, the Hankel functions have the asymptotic forms, 



Hi^^iz) ^ y^2/{^e-'^'~'^-^+^'^>'^\ (4.39) 



We set ci(fc) = y'7r/(4A;)e*'^'^"+^/^^'^/^, and C2(fc) = 0. The phase factor of ci(fc) is chosen in order to match with 
Eq. H4.29|l at short scales, while the factor of ^JiT/{'ik) ensures conformity with the Wronskian in Eq. H4.24|l . Therefore 
the final solution with the appropriate normalization is 

Quiz) = ^^^JWy^'"^^"^^^'^V~zHl^^{z). (4.40) 

It is worth mentioning that that the matrix valued Hankel functions are to be interpreted as series expansions, just 
like the usual Hankel functions. 

We finally discuss the solution in the super- horizon region. On super-horizon scales we have \z\ <ti 1, for which the 
asymptotic form of the Hankel function is 

' V / r(3/2) ^ ' 

so that the asymptotic solution for Qk{z) in the super- horizon region is given by 

Qu{z) ^ (i/V2fc)e*(^«-i/2)V22-K-3/2lV7^^p-.K^ 

i (o/2j 



-(l/V2fc)e*(^"-l/2+25„)V2^^(^/^^)-l-5K^ 



(4.42) 



where 

En = {l~€n)t + {2--iE~\n2)5n, (4.43) 

and ~ 0.5772 is the Euler constant. 

In this region since fc/7i ^ 0, we can also solve Eq. (I4.2t)|) ignoring the fc^ dependent term, leading to 

/"'' d-q' 

Ukiv) ^up k + CkO + Dk" 



6l2(ry')' 



Up k = 



-f / din"en7i^q2k, (4.44) 



where C'k and Dk are constants of integration, and up fc is a particular solution. Note that since is a rapidly decaying 
function, we can ignore Ck compared to Dk- In the same approximation, the solution of Eq. (j4.27|l is 

qik = dk{l/0) + 2(1/0) r dr]'enT]^q2k. (4.45) 

Jtj-h 

From Eq. 14.12|l we see that the integration constants Dk and dk are related by Dk = ^dk- Considering the region where 
T] is sufficiently close to rju, the integral in Eq. H4.45|l may then be neglected, so that using Eq. H4.32|l . we can write 
qik — 2D k{l / 0-}i){z / zt-c)~^ ■ Taking into account the asymptotic solution (|4.42l) . and the fact that qk = (ei • em)'^qik, 
we finally obtain, 

Dk = -(l/2V2fc)e'(''«-i/2+2««)^/20„(ei • e^fEn- (4.46) 

Thus the integration constant in Eq. (|4.44|) is completely determined to first order in slow-roll. Inserting the result 
for Dk in H4.6|l together with (|4.44|) . and using the relation a-uHu = k, we finally arrive at 

= e'(--V2+2.„)V2_^^ [yitnMei ■ e„,f + / {tn,t)] En, (4.47) 
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where we ignored Ck, and 

^{tn,t) = 2 f dt'H^^.[^^{e,-e^fR-^. (4.48) 
Jtn yea Q^n 

Here Qkn is the value of the asymptotic solution H4.42|l for Qk evaluated at rj — rj-H. Observe that the solution (|4.47() 
for (fi^k is expressed entirely in terms of background quantities and comoving time. 



4.5. Vector and tensor perturbations 



For the sake of completeness, we now present a brief discussion of vector- and tensor-type perturbations. From the 
G° component of Eq. H2.7(l . together with Eq. H3.13|l . we have 



(4.49) 



while the condition Tj^^ — yields 



2 



(4.50) 



Equations (|4.49|l and H4.5()|l describe the vector-type, or rotational perturbations. Observe that v^'"\ and tt*^^-' 

appearing in these equations are gauge-invariant, see Eq. (|3.18|) . Since vector sources are absent when the matter 
sector is composed entirely of scalar fields, the vector-type perturbations are therefore irrelevant in the inflationary 
scenario. 

The equation for the tensor-type, or gravitational wave perturbations follows from the G'^- component of H2.7|) : 



(4.51) 



For scalar fields we have tt^*^ = 0. Here Cl^ is symmetric, transverse-traceless, and gauge-invariant. We quantize the 
tensor perturbations by interpreting C*-*-* as the operator C**-*-* with the mode expansion 



(t) 



C, 



(fik 
(27r)3/2 



h.C. 



V — ey(fc; A)wfc^A(f7)- 



(4.52) 
(4.53) 



The quantity ey (fc;A) is the polarization tensor satisfying the conditions 



0, fc'e^j = 0, and 



e;.(fc;A)er(fc;A') 



The summation is over the two independent polarization states -I- and x . As before. 



the creation and annihilation operators satisfy the commutation relation. 



(4.54) 



The mode Vk^\{r]) satisfies 



a 
a 



Vk = 0, 



(4.55) 



for each A, and a" /a = 'H^{2 — e). Proceeding similarly as in the case of scalar perturbations, the solution of Eq. (|4.55|) 
in a sufficiently small interval of time around rj-j-c to the first order in slow-roll may be written in terms of Hankel 
functions (z = krj): 



(4.56) 



15 



Performing the asymptotic expansion of Hankel function in the region \z\ <C 1, the super-horizon solution for Vk[z) 
becomes 

Vk{z) ^ (l/A/2fc)e^(2e„+i)V2c^(^/^^)-i-.„^ 

Ch = 1 + (1 - 7£ - ln2)e„. (4.57) 
On the other hand, solving (|4.55(l in the super-horizon region yields 

1 

Vkiz) = Aka + Bka -^dz'. (4.58) 

We ignore the rapidly decaying term B^, and matching the term with the solution H4.57|l . we finally obtain the 
solution of tensor perturbation in the super-horizon region: 

clU^)-e^^'^^+'^''^'^CnHn E '^^i^^'^) (4-59) 

A= + ,x 

5. POST-INFLATIONARY PERTURBATIONS 

5.1. Adiabatic and Isocurvature perturbations 

In the previous Section we presented the super-horizon solutions to scalar perturbations during inflation. We now 
extend the analysis to density perturbations after inflation, in the regime of radiation- / matter-domination and 
recombination. 

It is well known that the most general density perturbation is a linear combination of adiabatic and isocurvature, 
or entropy perturbations. Adiabatic perturbations are just the total energy density perturbations. They perturb 
the solution along the same trajectory in phase-space as the background solution. On the other hand, isocurvature 
perturbations refer to the condition when there is no perturbation in the total energy density, there are perturbations 
in the ratios of energy densities of two or more of the components adding up to zero. Isocurvature perturbations 
perturb the solution off the background solution. When there is just one scalar field driving the inflation, isocurvature 
perturbations vanish identically. However, in the multicomponent situation, not only isocurvature perturbations are 
generated, but also they may be correlated with the adiabatic perturbations . According to the definition (|2.12|l 
of our basis vectors {e„}, when there are N scalar fields, the adiabatic perurbation will be along the vector ei, while 
the — 1 isocrvature perturbations will be along directions orthonormal to ei |19| |. 

In order to simplify our analysis, we assume that, of the N scalar fields, one of them has decayed to Standard 
Model particles so that isocurvature perturbations among them are absent, while the remaining A^ — 1 components 
have decayed to non- interacting Cold Dark Matter (CDM) js^. Then the CDM components may be considered as 
ideal fluids without mutual interactions. 

We start with Eq. (|3.34|) in conformal time for K = — tt^^\ with the gauge-invariant entropy perturbation e 
defined in Eq. (|XT7|) : 

^'^ + 3H{l + cl)^'^+clk\^ + [2H' + {l + 3cl)H']v^ 

= 2Kla^^icl - p)S, (5.1) 



where 



is the total entropy perturbation, and 



1 e _ Y^xLi^ t^jS j 



^ (5.3) 

Ml +Pi + Pi 

In Eq. 1)5. 3|l we used the notation that the Standard Model particles are represented as photons (with subscript 7) 
and Pi — for CDM components, while p^ = n-y/i for photons. Since we are considering the CDM components as 
mutually non-interacting ideal fluids, it can be shown that in flat space Si are constant in the super-Hubble region j37j . 
Further, /i^ have the same time-dependence. Therefore, within our approximation, S is constant for super-horizon 
scales. 
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5.2. Super-Horizon solutions 



The solution of Eq. H5.1|l in the super- horizon region, where we neglect the c^fc^ term, may be obtained in the 
same way that we derived the inflationary solution. Introducing the variables u and 9 defined in Eqs. (|4.6|) and (|4.11|) 
respectively, Eq. (|5.1|) becomes, 

u" - -rw = -2V3 ^='9, (5.4) 

Kq y/l+W 



with w — p/fi. This eqution admits the solution (p-^ — ip^ + f x j where 

^ _^2(^__3£)_ j dt'a{l + w) 



a 



to 



(5.5) 



is the homogeneous solution (with 5 = 0), and 



(^(^) = 6- f dt'a{l + w) f dt"H^^S 



' J 

(5.6) 

is the particular solution. The function is defined in Eq. (|4.48() . It immediately follows from Eq. (|5.6() that when 
(P) 

S is constant, (fx = 2S is a particular solution. 

Returning back to our discussion in Section l5. II we can claim that the adiabatic perturbation corresponds to the 
homogeneous solution of Eq. H5.1|l with the source term S absent, while the isocurvature perturbation pertains to 
the particular solution due to the entropic source S. The initial conditions that one imposes on the isocurvature 
perturbation are that both ip-^ and t/?^ vanish at the beginning of the radiation-dominated era. If we further make the 
simplifying assumption that the end of inflation at time te marks the immediate beginning of the radiation epoch, we 
can then write the solution to isocurvature perturbation with a constant S as 

^(^.o)(^) = 2s(l~2^^^^^^^^{t,,t) 



^ V l + w{te) 

= 25(^1 -^^(te,t)). (5.7) 

In the radiation-dominated era we have a{t) cx t^/^ and w = 1/3, so that the non-decaying part of J^{te,t) — 2/3. 
Thus, from Eq. (|5.7|l . t/^x*"^ vanishes for radiation domination. On the other hand, at the time of recombination 
during matter-domination, we have a{t) oc t^/^^ w = 0, and hence J^{te,t) = 3/5. In this case we have the simple 
relation 

f^r^ = (5.8) 

We now derive the expression for S in the super-horizon region during inflation. Using the slow-roll functions 
introduced in Section ITTl we have 

Z Kg 



w=- = -l + -e, c^ = - = -l--ry, 
/i 3 /i 3 

Sp-cl6fi^'^H^Vev^^ (5.9) 
Ko a 
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Substituting these into the definition S = \{5p — c^S^)/{p — c^/i) yields 

S^J^^r^^^. (5.10) 

From Eq. (|5.1()(l we see that the total entropy perturbation is generated along the directions orthonormal to the 
adiabatic perturbation, in agreement with Il3|. It is worth mentioning here that when e = 2He{e + 77") = 0, the 
expression for S as given by (|5.1U|) develops a singularity. Hence S is not a convenient variable to use during inflation. 
However, since we are demanding that the transition from the end of inflation to the radiation era is immediate, we 
will evaluate S at the time t,,- 

It remains to obtain the solution for the adiabatic perturbation. This is derived by matching the solution for 
ip-^ in Eq. (|4.47|) at the end of inflation t,,, with the homogeneous solution H5.5|) while maintaining continuity and 
differentiability. Ignoring the rapidly decaying Ck term, it follows that 



(5.11) 



The functions J and are defined in Eq. H4.48|l . 



5.3. Power spectra and spectral indices from inflation 

Having obtained the solutions for adiabatic and isocurvature perturbations given by Eqs. H5.11(l and (|5.8(l respec- 
tively, and the solution for the tensor perturbation H4.59I) . we are now in a position to calculate the power spectra 
and spectral indices from inflation. The power spectra are conventionally defined as 



AUk) ^ ^ (d^if 47^1 + l^rVff^l)) . (5.12) 

The subscripts A, S, and T refer to the adiabatic, isocurvature, and tensor power spectra, while A^(fc) is the spectrum 
due to cross-correlation between adiabatic and isocurvature perturbations. Substituting the results for <P^'fc'\ V'xfc"'' 
and Clj\ in the above expressions, and the using the fact that y{te,t) = 3/5 during matter-domination, we obtain. 



^Aik) = 4 fr!^ f(l - 2e„)(l + ^e'^^-) - 2C^o f(2ew + vli) + 2?7w(e2 • e^f^^ + X^Sn'^e)] , (5.13) 
A5(fc) = [(1 - 2e„)r/ye - 2C,y^5n%] , (5.14) 

A^(fc) = ^^^^^ [(1 - 2eH)^e'^^e - 2Co (77^(62 • emf% + ^1 5n^e)] , (5.15) 



AT(fc) = 16^ [1 - 2(Co + l)eH)] , (5.16) 



H; 

where Mpi = G^^l"^ is the Planck Mass, and 



^ = ^V^^^7T-(e2-e,„)^, Co = 7iJ + ln2-2. (5.17) 
2 e -t- 77II Qu 



The subscript e reminds us that and are evaluated at t — 
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If we now assume that the power spectra A^(fc) depend weakly on k, where X denotes adiabatic, isocurvature or 
cross-correlated spectra, we can parametrize them as 



A^(fc) = A^(fco)(^ 



rix(feo) — 1 



while the tensor spectrum is conventionally parametrized as 



AUk) = A|(fco) 



(5.18) 



(5.19) 



The normalization factor A^(fco) is called the amplitude, and n the spectral index. Here fco is a suitable pivot 
wavenumber. These parametrizations are valid for a range of k when nx — 1 and tit are close to zero, that is, when 
the power spectra are near scale-invariant. They lead to the definition of the spectral indices. 



dlnA|(fc) 
d Infc 



for the scalar modes and 



nrika) 



rflnAf,(fc) 
d Ink 



(5.20) 



(5.21) 



for the tensor modes, with the right hand side of Eqs. H5.20|l and 15.211 are to be evaluated a,t k — ko. Substituting the 
values of A^(fc) from above leads to the expressions for the spectral indices valid to first order in slow- roll: 



"-^ - 1 = - 

715-1 = -2 

nc-l = -2 

ut — — 2e-^ 



4ew - Zrjy^ H 



1 + 



(5.22) 

(5.23) 

(5.24) 
(5.25) 



The contribution due to multicomponet scalar fields come from and Y. In the case of a single scalar field, both 
these terms vanish, and we recover the single-field results. Observe that in this case ns and nc are irrelevant. 

We would like to point out here that one often comes across the adiabatic power spectrum defined in terms of 
the curvature perturbation TZ = —ip-^ + {H/H){ip^ + Hip^) and denoted by A^ [l3> US l23l ■ This definition simply 
removes the time-dependent factor in Eq. (|5.5(l . so that A'^ — (9/25)A^, and ua = n-ji- 

Using the above results we obtain the consistency relation, or the tensor to scalar ratio js^, 



A|(fco) 



16eH(l + ^e^'24)" 



(5.26) 



The single-field result r : 
may be written as [20ll3E 



16e-H 



■iuT follows immediately from above, while in the case of two fields, Eq. H5.26|l 



-8nTsin A, cosA 



A2 



A2 A2 



(5.27) 



For single-field inflation, WMAP hmits the tensor to scalar ratio as r(fco = 0.002 Mpc"^) < 1.28 (95% CL) 
An independent estimation of ut and r would provide a decisive test for discriminating multiple-field inflation from 
single- field slow- roll models The results obtained in this Section may be applied towards identifying classes of 
inflation models differing by their observational signatures l4l| . They mayalso be used to obtain information about 
the slow-roll potential by an inverse analysis of the observational results |39| . 
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6. CONCLUSION 

The foregoing Sections contain a general framework for analyzing the dynamics of cosmological perturbations 
in multiple-field inflation using the gauge-ready approach. The model comprises arbitrary number of scalar fields 
induced with a general field metric coupled to Einstein gravity. The complete set of equations describing scalar 
perturbations were presented in the gauge-ready form as well as in terms of gauge-invariant variables. Solutions for 
density perturbations were derived to the first order in slow roll during inflation, and for adiabatic and isocurvature 
modes after inflation. Tensor perturbations were also discussed. Expressions of the power-spectra and spectral 
indices for adiabatic, isocurvature, cross-correlated and tensor modes were obtained within the first order slow roll 
approximation. These results are of direct relevance when comparing theoretical predictions of various inflation models 
with observations. It would be of interest to interface our approach with the CMBFAST or CAMB |4^] computer 
codes and compare with the WMAP results. 
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